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, Abstract. We investigate the internal structure of cold dark matter 

H ' halos using high resolution N-body simulations. As the mass and force 

, resolution are increased, halo density profiles become steeper, asymp- 

toting to a slope of ~ —4/3 in the central regions and may not have 
j> . converged to a unique result. At our highest resolution we have nearly 3 

Q\ ' million particles within the virial radius, -R20O) an d force softening that is 

l/"") ■ ~ 0.2%i?2oo- This resolution has also allowed us to resolve a large part of 

the overmerging problem - we find over 1000 surviving dark halos orbiting 
within a single cluster potential. These data have given us unprecedented 
[ — \ insights into the dynamics and structure of "halos within halos", allow- 

Q^ ■ ing us for the first time, to compare the distribution of dark matter with 

observations of galaxies in clusters. 
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c/3 , 1. Introduction 

ctf ■ 

I We follow the path of many previous researchers who have used numerical sim- 

V~j ■ ulations to investigate the formation and structure of dark matter halos (c.f. 

Zurek et al. (1986), Frenk et al. (1988), Dubinski & Carlberg (1991), Warren 

1 et al. (1992), Carlberg (1994), Crone et al. (1994), Summers et al. (1995) 

Navarro, Frenk & White (1996, NFW), Tormen et al. (1996). The systematic 
study by NFW showed that halos density profiles may follow a universal form, 
uniquely determined by their mass and virial radius; varying from r _1 in the 
central regions, smoothly rolling over to r -3 at the virial radii. These studies 
have a resolution such that halos typically contained ~ 10, 000 particles and 
force softening of 1% of the virial radii. (The higher resolution simulation by 
Dubinski & Carlberg (1991) showed the same steep inner profiles, but had open 
boundary conditions such that infall ceased at z ~ 1. ) 

The CDM model has power on all scales, allowing most of the mass to 
collapse into very small halos at early epochs that are not resolved by numerical 
simulations. If this material at higher densities can be resolved, then it may alter 
the evolutionary history of the final halo, perhaps leading to different density 
profiles. One of the key aims of this work is to verify the convergence to a unique 
profile. 
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Moore, Katz & Lake (1996) pointed out that the amount of substructure 
that could survive within virialised regions depended simply upon the numerical 
resolution and that the overmerging problem resulted primarily from large soft- 
ening lengths. Halos in high density environments would be tidally disrupted 
since they are artificially "soft" . One of the results of the higher mass and force 
resolution is that halos can be found surviving within high density regions and 
thus for the first time we can examine the dynamics and properties of the dark 
halos within clusters. 

2. The simulations 

A candidate cluster is initially identified from a large cosmological volume that 
has been simulated at lower resolution. The particles within the selected halo 
are traced back to the initial conditions to identify the region that will be re- 
simulated at higher resolution. The power spectrum is extrapolated down to 
smaller scales, matched at the boundaries such that both the power and waves 
of the new density field are identical in the region of overlap, then this region 
is populated with a new subset of less massive particles. Beyond the high res- 
olution region the mass resolution is decreased in a series of shells such that 
the external tidal field is modeled correctly in a cosmological context. The new 
initial conditions are integrated with high accuracy using the parallel treecode 
PKDGRAV, using variable timesteps and periodic boundary conditions. We 
systematically increase the force and mass resolution for two halos with virial 
radii of 2 Mpc and 3.4 Mpc. At our highest resolution our force softening is 5 
kpc and the particle mass is 8 x 10 8 M© within a cluster with i?200 = 3.4 Mpc 
and mass ~ 2 x 10 15 M Q (we assume H=50 km s^ 1 Mpc -1 throughout). 

3. Results 

As we increase the mass and force resolution we find that the inner slope of 
the density profiles become as steep as -1.4. We suspect that the steepening of 
the slope from -1.0 found at lower resolution, results because we are resolving 
many more halos collapsing at early epochs when the mean density of the Uni- 
verse was high. This material is more robust against tidal disruption and can 
sink deeper within the cluster potential. Alternatively, Evans & Collett (1997) 
demonstrated that an inner slope of -4/3 is a stable solution to the Fokker-Planck 
equation and collisional Boltzman equation. 

Halos on galaxy scales will collapse earlier than the cluster simulated here, 
thus at a fixed radius measured in terms of the virial radius, galaxy halos will 
be even denser than clusters. These steep inner profiles will lead to problems 
when trying to reproduce the properties of giant arcs in clusters, the slow rise 
of stellar velocity dispersions with radius in Cd galaxies and the rotation curves 
of disk dominated galaxies. 

The distribution of halos orbits are close to isotropic; the median value 
of apocentric to pericentric distances is 5:1, a value that does not vary with 
position within the cluster and is unbiased with respect to the orbits of the 
smooth particle background. Circular orbits are rare and more than 25% of all 
halos within the cluster will pass within 200 kpc = 0.1i?200- 
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Figure 1. The upper panel shows the projected distribution of sub- 
structure halos that lie within i?200 • The solid circles give the measured 
tidal radii of the halos and the dotted circles denote their theoretical 
virial radii calculated assuming they are isothermal spheres. Mass loss 
via tidal stripping and halo-halo collisions is clearly evident. The length 
unit in this plot is in units of the cluster virial radius. Note that the 
covering factor of dark matter substructure is of order unity. The lower 
panel shows the inner density profiles of the entire cluster halo as we 
increase the force and mass resolution. The curves start at the values 
of the softening parameter used. As a rule of thumb, we use spline 
softening lengths that are 1 / 50th of the mean inter-particle separation 
within the whole box. Comparing the halo profiles at different reso- 
lutions demonstrates that we can "believe" a given profile to a scale 
of about 0.5 times the mean inter-particle separation of those parti- 
cles within the virial radius. Note that this is a smaller scale than 
advocated by Splinter et al. (1997) who argue that many statistics are 
affected on scales smaller that the global mean interparticle seperation. 
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Dark halos will be tidally truncated to a value determined by the density 
of the cluster at their pericentric positions. The approximation of isothermal 
halo mass distributions orbiting within a deeper isothermal potential works very 

Well; i.e. Tfidal = f 'peri° 'halos l&dus- 

The mass bound to resolved dark matter halos is approximately 10% of the 
entire cluster mass and varies from 0% near the cluster center, to 20% at its virial 
radius. This latter value is roughly the expected value for the mass attached to 
halos above a circular velocity of 80 km s _1 . (As numerical resolution increases, 
we might expect this number to asymptote to unity.) Correspondingly, through 
tidal stripping, the sizes of halos vary with cluster-centric radius, an effect that 
will be readily observable using weak lensing of background galaxies. 

Overmerging within the central regions of dense halos leads to a final distri- 
bution of substructure that is always anti-biased with respect to the global mass 
distribution. We find no evidence for a velocity bias - halos move on the same 
orbits as the smooth particle background and have the same velocity dispersion. 

The density profiles of a sample of well resolved halos within the cluster and 
those in the cluster proximity, have significantly higher concentrations than those 
found in isolated environments. This is due to their earlier collapse redshifts 
rather than the internal response of the halos to mass loss and heating from 
tidal stripping. 

Most of the halos within the cluster and in the cluster proximity have density 
profiles that are well fit by NFW profiles. Halos that loose a great deal of mass 
through tidal stripping have outer density profiles as steep as p{r) oc r~ 4 (at 
w 30% of their virial radius), thus Hernquist profiles provide slightly better fits. 

Mergers between halos with mass ratios larger than 10:1, occur with a 
frequency of about 20% in the cluster vicinity over the past 5 Gyrs. Once 
within the virial radius mergers are very rare; not a single merger occurs for 
halos within 90% of r2oo of the cluster center. 
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